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O ; Abstract 



In this last of three articles on the optical absorption of electrons in a half-filled 



Peierls-distorted chain we address the dimerized extended Hubbard model in 

a 

i ■ the limit of a large on-site interaction U. When the Hubbard interaction is 

g 

o 

o 



large both compared to the band width W and the nearest-neighbor interac- 
tion V the charge dynamics is properly described by the Harris-Lange model. 



^ ■ This model can be exactly mapped onto a model of free spinless Fermions in 

parallel (Hubbard-)bands of width W which are eventually Peierls-split. 

To determine the coherent absorption features at low temperatures we 
design and employ the "no-recoil approximation" in which we assume that 
the momentum transfer to the spin degrees of freedom can only be Aq$ = 
or Aqs = ir/a during an optical excitation. We present explicit analytical 
results for the optical absorption in the presence of a lattice dimerization 8 
and a nearest-neighbor interaction V for the Neel and dimer state. We find 
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that the coherent part of the optical absorption for V = is given by a single 
peak at u = U and broad but weak absorption bands for WS < \u> — U\ < W. 
The central peak at u = U only vanishes for 5 = in the Neel state. For an 
appreciable nearest-neighbor interaction V > W/2 almost all spectral weight 
is transferred to the Aqc = 0-exciton and the Age = 7r/a-exciton whose 
relative spectral weights very sensitively depend on both the lattice and the 
spin dimerization of the ground state. 



PACS1996: 71.10.Fd, 71.20.Rv, 36.20.Kd 
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I. INTRODUCTION 



In two previous articles (Gebhard, Bott, Scheidler, Thomas, and Koch I and II 1996) we 
have studied the optical absorption of non-interacting and strongly interacting tight-binding 
electrons in Peierls-distorted chains. These two different cases are idealizations of real al- 
most ideal one-dimensional semiconductors (Farges 1994). While polymers like polyacetylen 
(Heeger, Kivelson, Schrieffer, and Su 1988); (Baeriswyl, Campbell, and Mazumdar 1992); 
(Schott and Nechtschein 1994) are commonly believed to be Peierls insulators some charge- 
transfer salts (Alcacer 1994); (Brau and Farges 1994) should be understood as Mott-Hubbard 
insulators. The extended Hubbard model for interacting electrons on a distorted chain at 
half-filling is considered appropriate for the latter class of materials (Mazumdar and Dixit 
1986); (Fritsch and Ducasse 1991); (Mila 1995). Unfortunately, the study of the optical, i.e., 
finite frequency properties of correlated electron systems poses a very difficult many-body 
problem that cannot be solved analytically without further approximations on the Hamil- 
tonian or the calculations (Kohn 1964); (Maldague 1977); (Lyo and Galinar 1977); (Lyo 
1978); (Galinar 1979); (Campbell, Gammel, and Loh 1989); (Mahan 1990); (Shastry and 
Sutherland 1990); (Stafford, Millis, and Shastry 1991); (Fye, Martins, Scalapino, Wagner, 
and Hanke 1992); (Stafford and Millis 1993). 

In our second article (Gebhard et al. II 1996) we focused on the optical absorption in 
the dimerized extended Harris-Lange model which is equivalent to the Hubbard model to 
lowest order in the strong-coupling expansion. We could exactly map the charge degrees 
of freedom onto two parallel (Hubbard-)bands for free spinless Fermions of band-width W. 
The eigenstates of the Harris-Lange model are highly spin-degenerate which allowed us to 
exactly calculate the optical absorption. The results apply to the "hot-spin" regime of the 
Hubbard model when the temperature is large compared to the spin exchange energy J ~ 
W 2 /(U — V). Real experiments are carried out at low temperatures for which the system is 
in an unique ground state with antiferromagnetic correlations. Unfortunately, this problem 
cannot be solved analytically. In this article we will employ the analogy to an ordinary 
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semiconductor (electrons and holes in a phonon bath) to design a "no-recoil approximation" 
for the "chargeons" in a "spinon bath". It will allow us to determine the coherent absorption 
features of the Hubbard model at large U/W. 

The paper is organized as follows. In section [Tl] we derive our approximate Hamiltonian 
from the Hubbard model. The charge dynamics will be governed by the Harris-Lange model 
while the ground state is determined by the Heisenberg model. Basic results of our articles 
(Gebhard et al. I and II 1996) for the optical absorption and the current operator are 



recalled in section |Hll In section [TV] we introduce the "no-recoil approximation" which 
allows us to obtain informations on the coherent part of the absorption spectrum at low 
temperatures. The corresponding results for the Peierls-distorted extended Harris-Lange 
model are presented in section [V|. Summary and conclusions close our presentation. 

II. MODEL HAMILTONIANS 

A. Charge degrees of freedom at strong coupling: the Harris-Lange model 

As shown in (Gebhard et al. II 1996) the Peierls-distorted extended Hubbard model 
(Hubbard 1963); (Efiler and Korepin 1994) can be mapped onto the Harris-Lange model 
in the limit of strong correlations. In standard notation of second quantization the latter 
model reads 

£dhn >e xt = fhw{s) + fuHB(5) + U f) + V v (la) 
T LHB = (-*) EC 1 + (-1)'*) (1 - n^ a ) (c+ a c l+1>a + c+ +1> A,a) (1 - ^+w) (lb) 

l,cr 

Tuhb = (-*) £(1 + (-l^K-a (cf CT q +1)CT + ci +lja c U(T ) nj+i,-* (lc) 

D = ^2hi t ^hi,i (Id) 
I 

V = J2(n l -l)(n l+1 -l) (le) 
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where h^ a = Ci a c l a is the local density of a-electrons and hi = hi^ + is the total local 
electron density. UD is the Hubbard interaction between electrons on the same site, VV 
is the nearest-neighbor interaction between charged objects like double occupancies and 
holes, and T L hb (Tuhb) describes the motion of holes (double occupancies) in the lower 
(upper) Hubbard band. The Harris-Lange model is equivalent to the Hubbard model to 
order £(£/[/) -1 , t(t/U)°, and t(V/U)°. Due to the Peierls distortion the electron transfer 
between two lattice sites is modulated by ±t5. We are only interested in the half-filled case 
where the number of electrons N equals the even number of lattice sites L. 



B. Band structure interpretation 

For V = and in the absence of a lattice distortion the exact eigenenergies of the Harris- 
Lange model are the same as those of independent spinless Fermions moving in two parallel 
bands separated by U. For the case of linear optical absorption we may thus work with the 
effective band structure Hamiltonian 

H^T d = E [(U + e(k))ht + e(k)h l k ] (2) 

\k\<n/a 

with e(k) = —2tacos(ka), and h k = u^u k , h l k = l k l k for our fermionic quasi-particles 
(chargeons) in the upper and lower Hubbard band. Their band width is W = 4t, the lattice 
spacing is a, and the /c-values of the first Brillouin zone are spaced by 5k = 2 / k/(Lo). 

The lattice distortion results in a Peierls splitting of the upper and lower Hubbard band. 
In this case the effective Hamiltonian for linear optical absorption becomes (Gebhard et al. 
II 1996) 

^band = £ 1 {U + TE{k))fl U T + TE (k)h l J . (3) 

|Jfc|<7r/(2o),T=±l L J 

with the dispersion relation 



E(k) = \jt{k) 2 + A(k) 2 (4a) 
e{k) = -2tcos(ka) (4b) 



A(k) = 2t5 sin(ka) (4c) 

where n% in l k T ) are the number operators for the quasi-particles for the upper (r = +) and 
lower (r = — ) Peierls subband in the upper (u) and lower (I) Hubbard band. 

C. Spin degrees of freedom at strong coupling: the Heisenberg model 

At half-filling the ground state of the Harris-Lange model is 2 L -fold degenerate. We are 
ultimately interested in the optical properties of the (Peierls-distorted extended) Hubbard 
model at strong correlations for which the ground state for large U/W is unique. Conse- 
quently, we have to go to the next order in the expansion of the Hubbard model in W/U 
to lift the above degeneracy. The Harris-Lange model is unsatisfactory in yet another as- 
pect: at half filling a finite lattice distortion cannot be sustained within the model because 
there is no gain in electronic energy. The electronic part of the ground state energy is zero, 
irrespective of the Peierls parameter 5. 

A consistent treatment of the expansion in W/U will also produce corrections to the 
Harris-Lange model (Harris and Lange 1967); (van Dongen 1994) which would render the 
problem intractable. Thus we argue that the degeneracy of the ground state will inevitable be 
lifted by a residual spin-spin interaction which may have its origin in the itinerant exchange 
in the Hubbard model (next term in the expansion in W/U) or in their direct exchange which 
is not taken into account in the extended Hubbard model where only the direct Coulomb 
terms were kept ("Zero Differential Overlap Approximation" (Kivelson, Su, Schrieffer, and 
Heeger 1987); (Wu, Sun, and Nasu 1987); (Baeriswyl, Horsch, and Maki 1988); (Gammel 
and Campbell 1988); (Kivelson, Su, Schrieffer, and Heeger 1988); (Campbell, Gammel, 
and Loh 1988); (Painelli and Girlando 1989); (Campbell, Gammel, and Loh 1990)). The 
additional "Hubbard-H^-terms" might very well be important in strongly-correlated narrow- 
band materials. 

Thus we assume that the charge dynamics is still governed by the Harris-Lange model. 
In particular, we assume that the chargeons do not scatter from the spinons. For the spin 
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dynamics we choose the (dimerized, anisotropic) Heisenberg Hamiltonian with an anisotropy 
into the ^-direction 



^Hcis = £(i + (-i)'fc) \j± (s?s? +1 + sfsf +1 ) + j z s?st +l - - 
i L 4 



(5) 



where Si is the vector operator for the spin at site /, = Sf+iSf = c^q , , Sf = Sf—iSf = 
CyCj -p and S* = (n^ — hi^)/2. J± > and J z > are two generally different antiferro- 
magnetic coupling constants, and < 5s < 1 determines the degree of spin dimerization. It 
is well-known that a one-dimensional spin-system that is coupled to the lattice degrees of 
freedom shows the spin-Peierls effect, 5s > 0. 

If the itinerant exchange was responsible for the antiferromagnetic coupling we would 
have J = Jj_ = J z = (1 + 5 2 )2t 2 /U and 5 5 = 25/(1 + 5 2 ). In general, however, we do not 
assume a simple relation between Jj_, J z and t, or <5s and 5. We state, however, that Jj_, 
J 2 are supposed to be small energy scales compared to the Coulomb energies U, V. Its 
influence on the exact optical excitation energies will thus be neglected in this work. This 
implies, for example, that we cannot distinguish between singlet and triplet excitons. 

III. OPTICAL ABSORPTION AND EFFECTIVE CURRENT OPERATOR 

A. Optical conductivity and optical absorption 

The dielectric function e(u) and the coefficient for the linear optical absorption a(u) are 
given by (Haug and Koch 1990) 

e(u,) = 1 + (6a) 
n b c 

where Re{. . .} denotes the real part and rib is the background refractive index. It is sup- 
posed to be frequency independent near a resonance. Hence, the real part of the optical 
conductivity directly gives the absorption spectrum of the system. 
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The standard result (Maldague 1977); (Mahan 1990) for the real part of the optical 
conductivity in terms of the current-current correlation function x( u ) * s 

Re{a(c)} = Im{ * ( " )} (7) 

oj 

x(u) = ^i£dte**mU-) (8) 

where Af± is the number of chains per unit area perpendicular to the chain direction. 

The current-current correlation function can be spectrally decomposed in terms of exact 
eigenstates of the system as 

r i 1 

(9) 



u + (E n - E )+i>y oj- (E n - E ) + i>y _ 
Here, |0) is the exact ground state (energy E ), \n) are exact excited states (energy E n ), and 
|(0|j|n)| 2 are the oscillator strengths for optical transitions between them. Although 7 = + 
is infinitesimal we may introduce 7 > as a phenomenological broadening of the resonances 
at uj = ±(E n — E ). The spectral decomposition of the real part of the optical conductivity 
reads 

Re{a(cu)} = ^ £ I <°b»| 2 $ (w - (E n - E )) - 5 (u + (E n - E ))} (10) 

which is positive for all uj. 

In the following we will always plot the dimensionless reduced optical conductivity 

uRejaju > 0)} 

aUuJ >0)= M^W ■ 

Furthermore we replace the energy conservation S(x) by the smeared function 

W = 1 2I 2^ ( 12 ) 

7T{X Z + 7^J 

to include effects of phonons and experimental resolution. 



B. Effective current operator 



In (Gebhard et al. II 1996) we showed that the current operator for an excitation from 
the filled lower Hubbard band to the empty upper Hubbard band can be written as 
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&,.+ = -(«ca) £ (l + ("I)'*) (l + (-1)^ 

[ni +1 ,- ff c^i i(y c Ij<T (1 - - ni^ a cl a c l+l a (1 - n m ,_„.)] (13) 

where r] = —\Ri + i — R[ — a\/a < is the relative change of lattice distances due to the Peierls 
distortion (Gebhard et al. I 1996). j 1 I ^ ter+ creates a neighboring pair of a double occupancy 
and a hole. 

As a central result of (Gebhard et al. II 1996) we found the band structure representation 
of the current operator as 

i!Ter,+ = £ ^ae(fc)u+ +5/ J M/2 i ? . (14a) 

fe|,|g|<7r/a 

for the translational invariant case, and 

«band, dim 



iinter,'+ ~ E | ^ eCte (^) ^t+q/J-k-q/2 ^k+q/2+TT / J k-q/2+n / a ^ q (l^b) 

|g|,|fc|<w/(2o) k 



A (A;) 
o 



ii + 1 — ii + J 

a k+q/2+-K/a i k~q/2 a k+q/2 i k-q/2+-K/a 



X 



q+w/a 



for a Peierls-distorted lattice. The diagonalized form of the interband current operator is 
given in (Gebhard et al. II 1996). 

Since the excitation energy only depends on the charge configuration the complicated 
spin problem could be hidden in the operators x q which are defined only in terms of their 
product, 

1 



(6, r])x g ,(6, n)= £ - £ e**-**\l + r]S + (-1)' (5 + i7))(l +V5+ (-l) r (<5 + v)) 



S: 



^ T2 Z 

„ n I 



(0\S[,... SU (Ui+i - iih+x) , 5?, • • • 5i_ 2 > (15) 

{S'l-21 ■ ■ ■ S'r) (ir+ltr — tr+llr) > ^r-li - ■ ■ |0) . 

In practice, q' = q or q' = q + ir/a. In (Gebhard et al. II 1996) we studied the "hot-spin" 
case where all 2 L spin configurations were equally possible ground states. Now we will set 
up the "no-recoil" approximation for the case of the (unique) ground state of the Heisenberg 
model, equation @. 
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IV. NO-RECOIL APPROXIMATION 
A. Basic ideas 

Although the eigenenergies and even the eigenstates of the Hubbard model at strong 
coupling nicely display charge-spin separation (Ogata and Shiba 1990); (Parola and Sorella 
1990) the two dynamical degrees of freedom are coupled again by the current operator 
that creates an optical excitation. The oscillator strength for such an excitation, however, 
factorizes into a charge and a spin part. For the charge part our band structure picture 
applies. The spin part is determined by a ground state correlation function for nearest- 
neighbor singlets at all lattice distances, see below. Hence, the calculation of the optical 
absorption spectrum remains a difficult many-particle problem even for a formally charge- 
spin separated system. 

If we want to make further progress we have to make assumptions on the behavior of the 
spin system. In (Gebhard et al. II 1996) we treated the case of a degenerate spin background 
("hot-spin case") where all spin configurations were equivalent. Now we are interested in 
the more realistic case where the optical excitation starts from the unique ground state of 
the Heisenberg model, equation (|5|). 

Our band structure interpretation suggests an analogy between our strongly correlated 
Mott-Hubbard system and the situation in an intrinsic semiconductor with a direct gap. 
The Hubbard bands for the charges correspond to conduction and valence bands, while the 
spin degrees of freedom play the role of phonons in a semiconductor or metal. In both cases 
the energy scales are well separated, the spinon (phonon) energy being much smaller than 
the chargeon (electron) energy. The spin spectrum is gaped in the presence of dimerization 
(5s 7^ 0) such that the spin excitations correspond to optical rather than acoustical phonons 
in that case. This appealing analogy is also the basic idea behind the tJ-model approach to 
high-Tc superconductivity (Zhang and Rice 1988); (Dagotto 1994). 

In electron-phonon systems one may often ignore phonon emission or absorption during 
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an excitation at low temperatures. This is an approximation even at zero temperatures since 
phonons can always be emitted even if phonon absorption is impossible. Effectively there is 
a finite probability ( "Debye- Waller factor") that there is no momentum transferred to the 
phonons during an excitation by photons. This constitutes the "no-recoil approximation" . 
It has many applications in solid state science (e.g., X-ray scattering, MoBbauer-effect), and 
it is also very successfully applied for the explanation of optical spectra of semiconductors 
(Haug and Koch 1990). Thus we are confident that a (modified) "no- recoil approximation" 
will also be valid for the strongly correlated Hubbard model where chargeons and spinons 
are energetically well separated. 



B. Formulation of the approximation 

To formulate the approximation for our case we will scrutinize equation (13) for the 



spin sector. Recall that the charge sector has already been mapped onto a standard band 
structure picture. We apply 

\S[, . . . S , / '_ 1 , (|z, ii+x - I/, Tm) ) . . . S' L _ 2 ) = (16) 

A-(l-2)A — (!-2)m/ (\ \ I f \ Ql Ql \ 

1 S 1 S' \ b l U2, 4-3 — 4-2, I 3) j *->2) • • • 

for both expectation values in equation flI5|) where l~s shifts all spins by one site and Ts' 
shifts all spins but those at sites 2 and 3. 
Furthermore, we use the identity 

2 E (\ ~ S2S3) \S 2 , S 3 ) {S 3 , S 2 \ Q - S 2 S 3 ) = I (Ta, Is - Ts, i 2 ))(Us, U - I2, Ts) I • (17) 

S2,S3 

Now that we have eliminated the restrictions on the intermediate spin sum we can carry it 
out and are left with a complicated ground state expectation value, 

x+x gl = 2(0\Z+ 3 (q) Q - S 2 S 3 ) Z 2t3 (q')\0) (18a) 
Z 2 -M) = ^e-^ a ft 2 % (l - 2) (l + (-I)'*) (l + (-l)S) • (18b) 
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Thus far the expressions are exact. Note that almost -but not quite- the double 

Fourier transform of the correlation function of nearest-neighbor singlets at site I and r. 

To obtain further insight into the problem we first analyze two special cases which can 
exactly be solved before we formulate the general approximation. 

1. Neel state 

The Neel state is the ground state of the anisotropic Heisenberg model, equation fl5|), for 
J± = (Ising model), 

L/2 

|AF) = J] c+ / T c+ ma (vacuum) . (19) 
i=i 

It displays a perfect double-periodic structure, i. e., if we shift all spins by two lattice sites 
we recover |AF). Hence, 

ft 2 % (l - 2) \AF) = \[{l + (-I)') + (l - (-I)') f s f s , 1 ] |AF) . (20) 

The wave function renormalization factors become 

Z 2>3 (g)|AF) = ^±^|AF(g)) (21a) 



|AF( ff )> = ((1 + 5){l + V ) + e^(l - *)(1 - n)T s T sl l ) |AF) . (21b) 
Thus we may write 

\ (S q , + <W0 (V,o + <W/«) (AF(g)|i - S 2 S 3 |AF(g')> • (22) 



X q X q i 



AF 2 

For later use we define 



7 AF 
^0 



(AF(0)|^-S 2 S 3 |AF(0)> (23a) 



Zf = (AF(n/a)\- 4 - S 2 S 3 |AF(7r/a)> (23b) 



ryAF 



(AF(0)|^-S 2 S 3 |AF(7r/a))+h.c. 



(23c) 
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such that Xq x 
Explicitly, 

z AF = 2(5 + v y 

Zf = 2{l + 5r 1 f ■ (24) 

Ztl = A(5 + V )(1 + 5 V ) 
The Neel state dominantly provides the momentum Aq$ = vr/a to the charge system during 
an optical excitation since Z^ F > Zff > Zq F for generic values for 5, rj. 

We note that for 5 = r\ = only Z^ F contributes. This implies that for a Neel state one 
of the two Hubbard bands can be thought as being shifted by it /a which results in vertical 
transitions between antiparallel bands as for non-interacting electrons in a Peierls-distorted 
lattice (Gebhard et al. I 1996). 

2. Dimer state 

The ground state of the fully dimerized Heisenberg model, equation (|) for J± = J z = J 
and 5s = 1, is the dimer state 

L/2 fl 

| DIM) = II y 2 ( a it^+U ~ 4,|4+i,t) I vacuum) . (25) 

Again, the dimer state is invariant under a shift of all spins by two lattice sites. Hence, 
equations (p0|) , fl2"T|) , (p2|), and (|^) analogously hold for the dimer state. In particular, 

Q - S 2 S 3 ) |DIM) = |DIM) (26a) 
Q - S 2 S 3 ) f^DIM) = -i|DIM) (26b) 

holds such that the Z-factors for the dimer state become 
Z DIM = [l + or/ + 3(o + r/)] 2 /4 

Zf M = [3(l + 5r7) + 5 + r,f /A ■ (27) 

Zm M = [1 + SV + 3(5 + V)} [3(1 + Srj) + 5 + r,} /2 
We see again that momentum transfer by Aq$ = vr/a dominates in the dimer state but that 
there is zero momentum transfer even for 5 = r\ = 0, in contrast to the Neel state. 

13 



AF 



yAF Icy ~ + f> 
^0 /+ X 7T/a X 7r/a 



AF 



Z^ F /2, and [x+ /a x + x 



X n/a 



A F 



3. General case 



The above examples are limiting cases for the ground state of the general Heisenberg 
model in equation (g). The correlation function for nearest-neighbor spin singlets is long- 
ranged ordered such that the spin system only allows for momentum transfers Aq$ = 0,n/a. 
It has recently been shown by (Talstra, Strong, and Anderson 1995) that this order persists 
in the standard Heisenberg model (J± = J z = J, 6s = 0) as "hidden" long-range order 
although their is neither dimer nor antiferromagnetic order in the ground state. 

These observations make us confident that the following "no-recoil" approximation will 
give the coherent features of the optical absorption spectrum for all values of J±/J z > and 
5s in the dimerized Heisenberg model. 



nf s - l \o) 



,T + (-lV\ /i _(_!)»' 

w D w ( 2 J + Wbw [ 2 ' s s ' 



|0) +rest. (28) 



Due to the "hidden" long-range order of spin singlets in the ground state we expect that we 
need two different "Debye- Waller factors" wdw, »dw as in the cases of the Neel and dimer 
state where they both were unity. In general, the square of their absolute value is smaller 
than one as can be checked from sum rules. It is clear that more elegant approximations 
than this can be designed, e.g., an /-dependent factor u>dw(0 could be introduced to mimic 
a finite correlation length for finite temperatures as in (Gebhard et al. II 1996). We will not 
follow the latter ideas here. 

The wave function renormalization factors become 

Z 2 ,M\°} = 6q,0+ 2 9,n/a \n<l)} +rest (29a) 

|*(<7)> = (^dw(1 +5)(1 + V ) + e iga W^(l - <J)(1 - 11)%%^) |0> • (29b) 
Thus we may write 

^'L = I K° + 5 ^l a ) + - S^sl*^)) • (30) 

In the presence of a Peierls distortion there is no recoil from the spin system since the recipro- 
cal lattice vector is given by Q = ir/a. Even if the lattice distortion was absent, however, the 
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translational symmetry is broken by the "hidden" long-range order in the nearest-neighbor 
singlet correlation function. The coherent part of the spin contribution can then be expressed 
with the help of Xq x, 

V. OPTICAL CONDUCTIVITY IN THE HUBBARD MODEL AT STRONG 

COUPLING 

In (Gebhard et al. II 1996) we calculated the optical absorption for the case of an 
incoherent spin background. In the present case we only have to replace the expressions 
for the operator products x^x q , by those for the "no-recoil approximation", equations fl28[) 
and ©. 



noh 



Zo/2, 



X ir/a X ir/a 



roll 



Z n /2, and 



X -w/a X 



X n Ju 



^Tr/a 



coh 



Z M /2. 



A. Band case: V = 



Since the ground state implicitly breaks the translational symmetry of our system even 
for 5 = 7] = we directly approach the case of non- vanishing 8,7]. 

For the dimerized case we have to diagonalize the current operator, equation ( |14b|) , in 
terms of the quasi-particle operators for the respective Peierls bands. The result for the 
optical conductivity can be written as (Gebhard et al. II 1996) 

thA/1 



Re{cx(u; > 0,6,7])} 



\Kr'{k,q)\ 5(u - E T>T >(k,q)) (31a) 



t,t'=±1 | 9 |,|fe|<7r/(2a) 

with the absorption energies between the respective Peierls subbands 



E T y(k, q) = U + r'E(k + q/2) - rE(k - q/2) . 

The transition matrix elements are given by X +t+ (k,q) = 
\ + -(k, q), and 



(31b) 



-X--(k,q), \- + {k,q) 



\ +t+ (k, q) = iea 
\ + -(k, q) = iea 



A(k) 

-e(k)(a + {3*_ + {3 + a*L)x+ + -^(a + a*_ - P + (31)x+ +7v/a 



(32a) 
(32b) 
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Here we used the short-hand notation a± = ctk±q/2 etc. for the mixing amplitudes in the 
standard Bogoliubov transformation (Gebhard et al. I and II 1996). 

For the coherent part of the optical absorption we only need their values at q — 0. In 



this case, a\- fil = -e(k)/E(k) and 2a k (3 k = -A(k)/E(k). We thus find |A+,±(fc, q)\' 



Icoh 



S q ,o\K,±( k )\ 2 with 
|A +J+ (A;)| 2 

|A+,-(A;)| 2 



ea 



Zn „ + ArTTTTTTTTT + oZ 



'E(k)< 



{ea) 2 fe(k)A(k)X 



5 2 Z + Z n - 8Z_ 



e(k)A(k) \ 
ul 5E(k) ) 



M 



(33a) 
(33b) 



2 V 5E(k) J 

Each of the two quantities above can be expressed in terms of an expectation value of the 
projection operator (1/4 — S2S3). Hence we can be sure that the result is positive or zero. 
For vanishing nearest-neighbor interaction we may directly use equation (|31~1) to arrive 

at 



Re{cr coh (w > 0,5,//)} 



7T. 



[Jfe|<7r/(2o) 



(34) 



+T E \KA k )\ 2 (^-U-2E{k)) + 5{u-U + 2E{k)) 

|fe|<7r/(2a) 



For later use we define the abbreviations 
Hi 

Ho 



X{u) 



16 

W 2 L 


E 

|fe|<7r/(2a) 


E(k) 2 


16 
W 2 L 


E 

|fe|<7r/(2a) 


A(k) 4 
S 4 E(k) 2 


16 


E 

|fc|<7r/(2a) 


(e(k)A(k) 


W 2 L 


\ 5E(k) 



1 + 1*1 



1 16 
nW 2 L 



E 



uj-U 



k\<n/(2a) 



ik)A(k) V 

6E(k) J (u - U) 2 - AE(k) 2 



(35a) 
(35b) 

(35c) 

(35d) 



For vanishing nearest-neighbor interaction we only need the imaginary part of the func- 
tion X(u) which can easily be calculated. The real part of the coherent optical conductivity 
finally becomes 
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Re{cr coh (u; > 0,5,7/)} = ■ 



(z #i + S 2 Z n H 2 + 5Z A/ # 3 ) 5(cu - U) (36) 



16aa; 

2(5% + Z„ - SZ M )y/[(u - Uf - (W5f\ [W> - (u - uy\ 



7iW 2 (l-5 2 ) 2 \uj-U\ 

in the regions u = U and W5 < \uo — U\ < W . Thus the generic coherent features are (i) a 
5-peak at u — U which results from vertical transitions between the parallel Hubbard bands 
(Age* = Aq$ = 0), and (ii) a band-to-band transition feature for W5 < \u — U\ < W which 
results from vertical transitions between two antiparallel bands (Aqc = Aq$ = n/a). 

To check the sum rule we integrate over the reduced coherent optical conductivity and 
obtain 

/ ^d red>coh (^, 8, rj) = (Z + Z n ) (37) 
Jo W 16 

k DW | 2 (i + ^) 2 (i + v) 2 cr n + I^dwI 2 (i - s) 2 (i - v ) 2 cf 



where £7^ ven ' odd j s nearest-neighbor spin-correlation function on even/odd sites, see equa- 



tion ( |53| ) of the appendix. The comparison with the sum rule, equation (55), shows 
that |u>dw| 2 < 1 and |u>dw| 2 < 1 have to hold which indeed allows to interpret them as 
Debye- Waller factors. 

For 5 = 1, r] = (L/2 independent two-site systems) we have Z = Z^ = Zm/2 = 
8\wtw 1 2 C| ven . The contribution of the band vanishes again, and we obtain a single peak at 
uj = U with the coherent fraction |wdw| 2 = 1 of the total oscillator strength. 

The case of a vanishing Peierls distortion (8 = rj — 0) is particularly interesting. Recall 
that in the case of an incoherent spin background ("hot-spin" case) (Gebhard et al. II) we 
observed a logarithmic divergence at u — U. Now we find from equation (^) that 



Re{cr red ,coh(w > 0)} = — 

16 



WirZ 8(u - U) + 2Z v \l\ 



u-U 



(38) 



W 

in which the proper Z-values for 5 = rj = have to be inserted. In general, the absorption 
consists of a 5-peak at u — U and a semielliptic contribution which come from vertical 
transitions between the parallel (Aqc = Aq$ = 0) and antiparallel (Aq$ = Aqc = ir/a) 
bands. 
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For the Neel state Z^ F (5 = rj = 0) = 0, and Z^ F (5 = r] = 0) = 2, and the result of (Lyo 
and Galinar 1977); (Lyo 1978) is recovered. It is only the Neel state that suppresses the 
physics of the parallel Hubbard bands for 5 = rj = 0. If we had only considered this state 
as reference state, the physics of the parallel Hubbard bands would have been missed. The 
dimer state with Z® 1M (5 = rj = 0) = 1/4, and Z^ lM {5 = rj = 0) = 9/4 does show the generic 
features. The exact results for both states are depicted in figure |l|. 

For nonzero lattice dimerization (5 ^ 0, rj ^ 0) even the Neel state reproduces the generic 
situation. The coherent peak at u = U is not at all a consequence of the lattice distortion 
but a consequence of the parallel Hubbard bands. This has been overlooked in previous 
analytical investigations (Lyo and Galinar 1977); (Lyo 1978); (Galinar 1979). Numerical 
calculations (Campbell, Gammel, and Loh 1989) for 5 = rj = essentially give the results 
from the Neel ground state since this state dominates for small system sizes. Small traces 
of a peak at uj = U may have been washed out by the adopted smoothing procedure. For 
non-zero lattice distortion, however, the ^-resonance at u — U becomes clearly visible even 
in the numerical simulations. 

The exact results for the optical absorption of the Neel and dimer state in the presence 
of a lattice distortion are shown in figure 0. For realistic states the Z-factors will not be 
too different such that the figure should reproduce the generic situation for the Peierls- 
distorted Hubbard model at strong correlations and low temperatures. This observation 
further supports our "no-recoil approximation" . Figure has to be compared to figure 4 of 
(Gebhard et al. II 1996). It is seen that the linear absorption at the threshold u = U—W now 
shows a square- root behavior. This is also the case in the "hot-spin" case for large enough 5 
when the Peierls gap has not been smeared out. It is again seen that the significant features 
of the absorption spectrum have not changed much when we go from the Harris-Lange model 
to the strongly correlated Hubbard model. 



18 



B. Exciton case: V ^ 



In the presence of a nearest-neighbor interaction we have to solve an integral equation 
when the lattice is distorted. Note that only the absolute position of the resonances will 
be determined by the Coulomb parameter U/t while their relative position depends on V/t. 
For illustrative purposes we tune V/t independently of U/t and put aside the question of 
the stability of the ground state against the formation of a charge density wave. 

The calculations are lengthy and will not be redone here, see appendix B of (Gebhard 
et al. II 1996) for details. The result can be expressed in terms of two operator-valued 
functions Gi^q), 



dig) = itea fx+Fi(g) - x^ +n/a F 3 (q) 



G 2 {q) = itea -x+F 3 (q) + x^ /a F 2 (q) 



(39a) 
(39b) 



and three g-dependent functions ^1,2,3 (9) as 



Re{a(u > 0, V, 6,7])} = Re{a(u > 0, 5, 77)} + 



2WV, 



auo 



Im { E (i 



+ VF 1 )(1 + VF 2 ) 



{VF,Y 



G x G^ 



G 2 G 2 



+ V (G x G+F 2 + G 2 G+F, + (G 1 G+ + G 2 G+)F 2 



(40) 



In the "no-recoil approximation" we only need the q = value of the functions -^1,2,3- We 
have 



F 1 (q = 0) 
F 2 (q = 0) 
F 3 (q = 0) 



#1 



uo-U 
h 



S 2 H, 



uj-U 



LO-U 



tt5 2 X(uj) 
ttX(lu) 
+ tiX(lu) 



(41a) 
(41b) 
(41c) 



see eqs. 



The definitions ( P3|) allow us to reduce equation ( f40|) to 
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Re{a coh (u > 0, V, S, rj)} = Re{a coh (u > 0, 5, 77)} 

VN±_{Wea> 2 ( 



16au lm \(l + VF 1 )(l + VF 2 )-(VF 3 y (42) 
Zo(F 2 + F 2 ) + Z W (F| + F 3 2 ) - Zm(F\ + F 2 )F 3 

+ V (F 1 F 2 - F 2 ) {Z F 1 + Z^F 2 - Z M F 3 



Note that the positions of the excitons are determined by the zeros of the denominator 
which is a function of u, V, and 5 but is independent of the Z-factors and rj. The oscillator 
strength of an exciton, however, strongly depends on the Z-factors. 

For the translational invariant case, 5 = 77 = 0, equation ( fl2|) can be simplified to (_F 3 = 0, 
F x = l/{u-U), F 2 = -irX{u)) 

R/r sni/n (Wea) 2 Af ± J \ ZpF, Z n F 2 \ 

From the first part we see that the peak at u = U for V = becomes a q = 0-exciton 
at frequency u> = U — V . For the Neel state this exciton is absent while it is present in 
all generic cases, e. g., in the dimer state. The second part describes the formation of the 
q = 7r/a-exciton from the semielliptic "band absorption" part for V — 0. The formation of 
this exciton leads to a redshift in the "band absorption" part as already noticed by (Galinar 
1979). 

The exact result for 5 = 77 = for the dimer state is displayed in figure ^| for various 
values of V/t. For V > W/2 the q = 7r/a-exciton becomes a true bound state at u w / a = 
U — V — W 2 / (4V) below the band absorption edge at u; e d ge = U — W. For large enough V/t 
the fact that Z n > Zq implies that the q = 7r/a-exciton resonance is stronger than the one 
from the q = 0-exciton. Since they are only separated by an energy difference 5uo = W 2 / (41/) 
thermal and disorder broadening could merge the two exciton lines into a single asymmetric 
absorption line. 

The general case of the optical absorption of the dimer state in the presence of a lattice 
distortion is shown in figure f|. The Neel state results in a very similar curve such that 
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we are confident that it reproduces the generic features for the coherent absorption in the 
Hubbard model at strong correlations. At V = we had the resonance at u — U and 
a Peierls-split "band absorption" in the range W5 < \cu — U\ < W. For moderate V/t 
the two bands evolve into two peaks such that three lines are visible. The peak at energy 
^ir/a = U — V — W 2 /(4V) corresponds to the q = 7r/a-exciton while the dominant peak is 
the q = 0-exciton at ujq = U — V . The peak near u = U corresponds to the anti-bound 
exciton which quickly looses oscillator strength as V/t increases. For large V/t the q = 0- 
and q = 7r/a-excitons dominate the absorption spectrum with a preference of the q = ix/a- 
exciton since Z n > Zq. In addition, a weak structure remains near uj = U which does 
not play a role for the optical absorption but may become visible in electroabsorption, i.e., 
optical absorption in the presence of a static electrical field. 

The dimer state underestimates the factors Z w , Zm as compared to Zq. For the Neel state 
the strength of the q = 0-exciton is much weaker for realistic values for the lattice distortion. 
The result for the Neel state is shown in figure [5[ It is now seen that the q = 7r/a-exciton 
strongly dominates the q = 0-exciton. This is to be expected since the components which 
imply a momentum transfer q = to the spin system (Zq, Zm) are much weaker now. Even 
a large value of 6 does not drastically change this situation. It should be clear that the 
strengths of all exciton peaks above the exciton at oo^/a = U — V — W 2 / (4V) are a direct 
measure for the dimerization degree both of the charge (5, rj ^ 0) and the spin (5s ^ 0) 
system. 

VI. SUMMARY AND CONCLUSIONS 

In this paper we have further studied the optical absorption for strongly-correlated elec- 
trons in half-filled Peierls-distorted chains. In (Gebhard et al. II 1996) we exactly solved the 
case of the Harris-Lange model with its incoherent spin background ("hot-spin case"). For 
the more realistic situation of low temperatures, i. e., a unique ground state, exact solutions 
are impossible for the generic cases. To make further progress we assumed that the spin 
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dynamics is irrelevant: (i) we neglected corrections of the order J to the optical excitation 
energies and thus a singlet/triplet-splitting of the excitations, and (ii) we assumed that no 
spin excitations are created during an optical excitation ("no-recoil approximation"). Then 
we could analytically investigate the problem even in the presence of a lattice dimerization 5 
and a nearest-neighbor interaction V between the electrons. 

In the "no-recoil approximation" for the strongly correlated Hubbard model all the fea- 
tures already seen in the "hot-spin case" become more prominent since now only transitions 
with Aqs = and Aqs = n/a are allowed. This corresponds to vertical transitions between 
parallel Hubbard bands (Aqs = Aqc = 0) and antiparallel bands (Aqs = Aqc = n/a). 
The relative strength of the two processes can be measured by three parameters (Z , Z n , 
Z M ) which depend on both the spin structure of the ground state and the lattice distortion. 
In the generic case the parallel Hubbard bands result in a single peak at uj — U . The 
Aqs = Aqc = ii /a part supplements the optical absorption spectrum by two absorption 
bands for WS < \uo — U\ < W . An additional nearest-neighbor interaction results in a 
tightly bound exciton at u = U — V and another exciton at u w / a = U — V — W 2 /(4:V). 
Simpson's exciton band has seized to exist since the exciton's center of mass momentum has 
to be Aq c = or Aq c — n/a. Again, the excitons draw the oscillator strength from the 
"band absorption" part. Although the location of the exciton resonances is independent of 
the spin configuration of the ground state their relative amplitude is a very sensitive mea- 
sure of the lattice and spin dimerization. We thus come to the important conclusion that 
the linear optical absorption in strongly correlated electron systems might serve as a subtle 
probe for the magnetic structure of the ground state. In particular, optical absorption allows 
to measure the "hidden" long-range order of nearest-neighbor spin singlets (Talstra, Strong, 
and Anderson 1995). 

One might think that it is fairly simple to experimentally distinguish between a Peierls 
and a Mott-Hubbard insulator on the basis of their optical absorption. At low temperatures 
we expect a Peierls insulator to show a broad band-to-band transition while a Mott-Hubbard 
insulator with appreciable values for the Coulomb interaction (U ^> W, V > W/2) should 
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display sharp exciton lines in the optical absorption spectrum. Unfortunately, a direct com- 
parison to experiment is difficult for two reasons. First, disorder effects can inhomogeneously 
broaden single lines. Hence, an experimentally observed band can very well be a sign of dis- 
order rather than an argument against a Mott-Hubbard insulator. Secondly, when a residual 
electron-electron interaction is included in a Peierls insulator (Abe, Yu, and Su 1992); (Abe, 
Schreiber, Su, and Yu 1992) one can equally well obtain excitons which draw the oscillator 
strength form the band transitions. Hence, single exciton lines are not a clear-cut indication 
against the presence of a Peierls insulator either. 

It is clear that one needs to compare model calculations and experimental observations 
on many more physical quantities than just the linear optical absorption before definite 
conclusions can be reached. 
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APPENDIX: SUM RULE 



We briefly account for the sum rule. We have 



r°° A/ 1 i ° \f\ 

/ du lm{ X {uj)} = k-±Y: W» =4^<0|j 2 |0>. 
Jo La „ 1 La 



It is a standard exercise to show that 



(0|j 2 |0) = (2tea) 2 £ (l + (-l)^) 2 (l + (-l)^) 2 (0| (± - SjSj+i ) |0> 



We define the positive quantities 



CT n ' odd = \ E^y^(°l Q - S.S Z+1 ) |0> 



(Al) 



(A2) 



(A3) 



such that Cs = Cg dd +Cg dd is the value of the nearest-neighbor spin-spin correlation function. 
Then we may write 

r du Im{xM} = nM ± a(2te) 2 \(1 + 5) 2 (1 + r]) 2 C7 cn + (1 - 5) 2 (1 - r,) 2 C° s dd ] . (A4) 
Jo 1 J 



The area under the curves for a ie d(oj), equation (pTT]) , is thus given by 
du 



t \ 11 
a rcc l ( u ) — — 

o W V ; 4 



(1 + 5) 2 (1 + t]) 2 C| ven + (1 - 5) 2 (1 - t]) 2 C° s dd 



(A5) 
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FIGURES 

FIG. 1. Reduced optical conductivity, a m d(o~> > 0), for the dimer and the Neel state for U = 2W. 
A broadening of 7 = 0.01IF has been included. 

FIG. 2. Reduced optical conductivity, <T re d(<^ > 0,5,r]), for the dimer and the Neel state for 
U = 2W in the presence of a lattice distortion, 5 = 0.2, rj = —0.06. A broadening of 7 = 0.01W 
has been included. 

FIG. 3. Reduced optical conductivity, <r re( j(a; > 0, V), for the dimer state for U = 2W in the 
presence of a nearest-neighbor interaction, V = 0, W/2, W. A broadening of 7 = 0.01IF has been 
included. 

FIG. 4. Reduced optical conductivity, o- rc d(<^> > 0,V,S,rf), for the dimer state for U = 2W 
in the presence of a lattice distortion, 5 = 0.2, i] = —0.06, and a nearest-neighbor interaction, 

V = 0, W/2, W. A broadening of 7 = 0.01W has been included. 

FIG. 5. Reduced optical conductivity, <7 re d(w > 0,V,6,rj), for the Neel state for U = 2W 
in the presence of a lattice distortion, 5 = 0.2, i] = —0.06, and a nearest-neighbor interaction, 

V = 0, W/2, W. A broadening of 7 = 0.01IF has been included. 
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